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Abstract 

The functional dependence of the high-energy observables of total cross sec- 
tion and slope parameter on the sizes of the colliding hadrons predicted by 
the model of the stochastic vacuum and the corresponding relations used in 
the geometric model of Povh and Hiifner are confronted with the experimen- 
tal data. The existence of a universal term in the expression for the slope, 
due purely to vacuum effects, independent of the energy and of the particular 
hadronic system, is investigated. 

PACS Numbers : 12.38 Lg , 13.85 Dz , 13.85 Lg . 
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1. The Model of the Stochastic Vacuum and the Geometric Models 

Diffractive high energy scattering is largely determined by the nonperturbative regime 
of QCD. The extended character of the interaction, involving correlation properties of the 
gauge field, determines the phenomenological properties of the observables, which are fixed 
by the sizes and global structures of the colliding systems, rather than by the number 
of their pointlike constituents and their couplings. These features have led to models of 
geometric nature for high-energy scattering which give natural account of the relations 
between total cross sections of different hadronic systems and, through hadronic form factors 
appropriately introduced, describe the shapes of the diffractive peaks. QCD must provide 
the fundamental framework in which this phenomenology should arise naturally, and efforts 
have been made in this direction. Some treatments based on perturbative QCD have also 
lead to dependence of the observables on hadronic sizes ||. 

A nonperturbative QCD description of the main features of high-energy scattering is 
given by the model of the stochastic vacuum |4]||, which combines QCD quantities (gluon 
condensate and corrrelation length) and hadronic sizes in an eikonal framework, leading 
to a unified description of the data for different hadronic systems. The calculations lead 
to definite size dependence of the observables of total cross section <j T and forward slope 
parameter B, and provide a dynamical scheme that explains the intuitive and Regge-based 
geometric treatments. In the present work we explore further the results of the model for 
the pp, pp and other hadronic systems, and compare them with those from the geometric 
model of B.Povh and J.Hiifner 0, confronting all predictions with the experimental data. 

In the calculation with the model of the stochastic vacuum the hadronic structures 
enter in the form of transverse wavefunctions (two-dimensional wavefunctions in the plane 
perpendicular to the direction of the colliding hadrons). Taking into account the results of 
the previous analysis of different hadronic systems |3J] we here only consider for the proton 
a diquark structure where, with respect to the relevant color degrees of freedom, the proton 
is described as a meson, in which the diquark replaces the antiquark. Thus we can treat on 
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equal footing meson-meson, meson-baryon and baryon-baryon scattering. Other structures 
for the baryons have been explored j|,|6| in investigations with the model of the stochastic 
vacuum, and the similarity of results obtained with different structures demonstrates the 
role of the extended nature of the nonperturbative dynamics. 

For the hadron transverse wavefunction we take the simple ansatz 

MR) = v^ir exp (~ r2 / s h) > C 1 ) 

where Sh is a parameter for the hadron size. The dimensionless scattering amplitude Th 1 h 2 
is given in terms of the dimensionless profile function Jh 1 h 2 f° r hadron-hadron scattering by 

T HlH2 = is[(g 2 FF)a A ] 2 a 2 J d 2 b exp (iq- b) J HlH2 (b, S u S 2 ) , (2) 

where the impact parameter vector b and the hadron sizes Si , 5*2 are here written in units of 
the correlation length a, and q is the momentum transfer projected on the transverse plane, 
in units of 1/a, so that the momentum transfer squared is t = — \q\ 2 /a 2 . For short, from 
now on we write Jib) or Jib/ a) to represent Jh^^P, Si, S 2 ). The normalization of T HlH2 is 
such that total and differential cross sections are given by 

a =-lmT HlH2 , — = —\T BlBi \ . (3) 

The observables are written in terms of dimensionless moments of the profile function 
(as before, with b in units of the correlation length a) 



d 2 bb k J(b) , k = 0,1, 2, ... (4) 
which depend only on Si/ a, S 2 /a, and the Fourier-Bessel transform 

I{t) = J d 2 b J Q (ba^\t\) J{b) , (5) 
where Jo(ba^J~\t\) is the zeroth-order Bessel function. Then 

T Hi h 2 =is[(g 2 FF)a 4 } 2 a 2 I(t) . (6) 
Since J(b) is real, a T and the slope parameter B are written 



d L da e£ \ 



= I h a 2 = Ka 2 . (7) 



It is important to observe that these results conveniently factorize the dimensionless 
QCD strength (g 2 FF)a 4 in the expressions for the observables. The correlation length a, 
which is an intrinsic parameter of the correlation function of the QCD field, appears as the 
natural length scale for the observables and for the geometric aspects of the interaction. 
These aspects are concentrated on the quantities Io(Si/a, S^/a) and ^{Sx/a^S^/a), which 
depend on the hadronic structures. These quantities are mainly determined by the values of 
the profile functions in the range of impact parameters up to about 2.5 fm. a T measures the 
strength, while the slope B has the strength cancelled out and is only related to the hadron 
geometry. The explicit formula for the slope is 

2 fd% J(b) 2 w v ; 

where it is seen as related to the average value of the square of the impact parameter in the 
collision, with J(b) as weigth function. We recall that here b is dimensionless and that (b 2 ) 
depends on the hadronic sizes. 

2. pp and pp systems 

We first discuss pp and pp systems, with Sx = 5*2 = S. The curves for I = 
a T / [{g 2 FF) 2 a 10 } and K = B/a 2 can be parametrized as simple powers of S/a with good 
accuracy, the convenient expressions being 

Jo = a(-) , K = v + 1 - . (9) 



a J ' \ a , 

The values of the parameters result from integrations over correlation functions |4[], are 
intrinsic to the model of the stochastic vacuum, and do not contain any dependence on 
experimental quantities. For the present purpose of analysis of data in a limited energy 
range, and for easier comparison to the geometric model we take their values as rj = 2.03, 
P = 8/3, 7 = 3/8, 5 = 2 , 6//3 = 3/4 , a = 0.76 x 10" 2 . 
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The proton radius can be eliminated from Eqs. fl7|) and (|]), and we obtain a relation 
between the observables a T and B at a given energy 

a 2 7 ( '<J T 

{B -^ )= [< g *FF>aT'^ V 2 ^) • (10) 

The two QCD parameters, (g 2 FF) and a, can be determined using this expression and the 
experimental data for o T and B at two different energies. 

The available data on a T and B in pp and pp scattering at high energies consist mainly 
of ISR (CERN) measurements at energies ranging from \/s = 23 GeV to y/s = 63 GeV, of 
the -v/i = 541 — 546 GeV measurements in CERN SPS and in Fermilab, and of the y/s = 1800 
GeV data from the E-710 Fermilab experiment. The Fermilab CDF measurements |§ at 
y/s = 1800 GeV seem discrepant with the E-710 experiment at the same energy and are not 
used here. 

Since we are here concerned with nonperturbative contributions only, at the ISR en- 
ergies we take for total cross sections the values given by the Donnachie and Landshoff 
parametrization || for the pomeron-exchange contribution 

< om (PP,PP) = (21-70 mb) S °- 0808 , (11) 

and for values of the slope we take those of the pp system (not those of the pp data). 
Using as input the data for the highest energies (541 and 1800 GeV), where the process is 
essentially nonperturbative and no separation is needed, we obtain || 

a = 0.32 ± 0.01 fm , < g 2 FF > a 4 = 18.7 ± 0.4 , < g 2 FF >= 2.7 ±0.1 GeV 4 . (12) 

The relation between the experimental values of the two observables is well represented 
at all energies from 23.5 to 1800 GeV with the form 

B = B A + C A (a T ) A . (13) 

We use this expression with B and Ba in GeV 2 , a T in mb and C A in mixed units. This form 
is similar to Eq. (|ToD , with an obvious correspondence of parameters. In our calculation 
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with the model of the stochastic vacuum the exponent A = 5/(3 does not depend on QCD 
quantities and is equal to about 0.75. 

Following ideas that relate hadron-hadron scattering to the shape and size of the colliding 
hadrons, B.Povh and J.Hiifner show that the combination of Regge amplitudes with 
electromagnetic form factors relates the slope parameter and the sum of the squares of the 
radii of the colliding hadrons, and write for Hp (hadron-proton) scattering 

B Up = Rl + K , (14) 

which is to be considered as a definition of effective hadronic radii. Observing the behavior 
of the experimental points in a plot of the observables a T and B against each other, they 
suggest that the dependence of the total cross section on the radii is 

a^ p = gR 2 R R 2 p . (15) 

Introducing electromagnetic form factors to reproduce the shape of the elastic differential 
cross section, B.Povh and J.Hiifner write the relation of these hadronic radii to the electro- 
magnetic radii as < r 2 m >= 3R 2 . Deviations from these simple formulae occurring at low 
energies show that they should be used only for y 7 ? > 20 GeV. 

For proton-proton scattering, with i?n = R P , Eqs. ([14] ) and ( |15|) lead to 

B = C 1/2 (O 172 , (16) 

which is of the form of Eq. fll3|) with A = 1/2 and B& = 0, and it is remarkable that, using 
as input the data at 541 and 1800 GeV, one obtains with A = 1/2 the same value B& = 0. 

In Fig. 1 we show the description of the data through Eq. (0), using for A the values 1, 
0.75 and 0.5. The case A = 1 is included in Fig. 1 for a numerical reference, although we do 
not refer to any model suggesting it. The vertical axis represents the constant Ca, which, 
together with Ba, is fixed in each case by the input data at the energies 541 and 1800 GeV. 
Then the five ISR data points are considered as parameter-free predictions, and we calculate 
a x 2 value representing the observed deviations. In the horizontal axis we mark the energy, 
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which here works just as an external label used to spread the information in the plot. There 
are no free parameters, since Ba and Ca are fixed by the input data. The values of x 2 
are also shown, and, although A = 0.75 is favoured, we cannot say that the differences are 
statistically meaningful. However, the model of the stochastic vacuum gives precise meaning 
to the parameters Ba and Ca in terms of QCD quantities, successfully predicts A = 0.75 , 
and introduces hadronic sizes in definite form, as parameters accounting for the extensions 
of the wavefunctions. It is remarkable the presence in this case of a bounding minimum Sa 
(equal to rja 2 ) for the slope, which is the same for all hadron-hadron systems. 



Fig. 1 - Test of the parameters of Eq. ( |l3| ) , comparing different models, using the experimental 
quantities of the pp and pp systems. For energies up to 62.3 GeV the values of a T are given by the 
parametrization 0"p Om = (21.70 mb)s 0808 and the values of B are those of the pp data. The values 
of i?A and Ca are obtained with the 541 GeV and 1800 GeV data as inputs. The horizontal lines 
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represent the constant Ca with the choices A=l, 0.75 (model of the stochastic vacuum) and 0.5 
(geometric model), x 2 represents the average deviation of the five ISR points from the constant 
line. 



The description given in Fig. 1 covers all pp and pp data. Extrapolating Eq. (|i3|) to 
higher energies (e.g. LHC energies), where the total cross sections may be about 100 mb, we 
find a small, but hopefully measurable, difference in the values of the slope, with B higher 
by 0.2 GeV -2 for A = 3/4 , compared with the A = 1/2 case. 

The proton radius presents a slow increase with the energy, taking values about the 
electromagnetic radius. In the case of the model of the stochastic vacuum, where the radius 
enters as a parameter of the wavefunction, the energy dependence of the radius can be 
parametrized in the form |J 

S p ( s ) = 0.671 + 0.057 log Vs (fm) (a), 
or S p (s) = 0.572 + 0.123 [log v^] - 75 (M (6). (17) 



In the geometric model Eq. (O) requires a power 1/2 in the logarithm in order to yield a 



log 2 s dependence in the cross section. 

Parametrizations (a) and (b) predict for ^/s = 14 TeV values of the proton radius 1.215 
fm and 1.240 fm respectively, which are about 40% higher than the electromagnetic radius. 
The pp cross section at this energy is predicted as (a) 95.5 mb and (b) 100.8 mb. These 
values are in good agreement with the results of the Akeno collaboration fl~2|j . 

The value of g is obtained by B.Povh and J.Hiifner, on the basis of Hp data, as g=75 
fm -2 . In Fig. 1 we show that C1/2 ~ 2 GeV~ 2 mb -1 ^ 2 , corresponding to g — (2/C1/2) 2 x 
0.1/(0.197) 4 = 66.4 fm" 2 . This means a fair simultaneous description of pp and pH systems, 
which is also achieved in the calculations with the model of the stochastic vacuum, where the 
interaction strength and the hadronic radii appear as fundamental quantities. An important 
phenomenological difference between the two approaches rests in the existence of a finite 
universal minimum value for the slope, which may be tested in systems where smaller 
hadrons collide with the proton, as we discuss below. 
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Eqs. (|14]), (p~5|) and (|16|) for pp and pp scattering can be built from a profile function of 
simple Gaussian shape 



j(b) = ^ e- b2 ' m l , (IS) 



_J,2 I a f?2 



where R p (s) has an energy dependence, and where the normalization for J is chosen appro- 
priately. However, two Gaussians are needed to describe the data for t ^ 0. 
3. Hadron-Proton Systems 

We now consider other systems of hadrons colliding at high energies. In the treatment 
of the pp system we are constrained by ^/s > 20 GeV, and cannot observe clearly the effect 
of the minimum slope . The contribution of this term could be better observed in Hp 
systems, where H represents hadrons of small size. We must remark that, since we deal with 
radii which are energy dependent quantities, we must compare different hadronic systems 
at the same center-of-mass energy. 

The parametrization of the results obtained with the model of the stochastic vacuum for 
general Hp systems is 

°Ln = h [{9 2 FF)a^ = [{g*FF)ara 2 a(^Y /2 , (19) 



a a 



and 



B = \j* 2 =ria 2 + \ 1 (sl + Sl) . (20) 



'0 

With a = 0.32 fm, we have r/a 2 = 5.38 GeV -2 . In the treatment of B.Povh and J.Hufner 
the corresponding relations are given by Eqs. (|14D and (|T5|) . 

In order to compare the models, it is important to eliminate the influence of specific values 
of radii, since they have different definitions. Thanks to the convenient factorization in the 
final expressions, we may actually build relations involving only the observables, or involving 
only the ratios of radii, which we may assume to follow the ratios of electromagnetic radii. 
We thus have for crHp/cr H 'p the ratios (th/th') 4 ^ 3 and (r H /rH') 2 m the stochastic vacuum 
and geometric models respectively. Entering with the known values [|l]]] for the radii of the 



proton (0.862 ± 0.012 fm), of the pion (0.66 ± 0.01 fm) and of the kaon (0.58 ± 0.04 fm) we 
obtain the results shown in Table I. The experimental ratio refers to the pomeron exhange 
contribution, taken from the parametrization of Donnachie and Landshoff. We observe that 
the value 2/3 given for the ratio <J n p/a P p by the quark additivity rule is here obtained as 
a simple consequence of the sizes of the hadrons. Also the ratio <7Kp/&irp is consistently 
obtained with a value close to the data, without need for different couplings of the pomeron 
to strange and non-strange quarks, as must be the case with quark additivity rules. Since 
the analysis of the proton structure in HERA (DESY) shows that the proton is better 
characterized as a sea rather than as a valence structure, the explanation of the high-energy 
phenomenology through the hadronic sizes is more legitimate. The factorization relation 
(j nn = a 2 p / (jpp is identically satisfied in both cases considered here. 

Table I - Ratios of the pomeron exchange contributions to total cross sections for different 
hadronic systems. The experimental values are taken from the parametrization of Donnachie and 
Landshoff. 



Cross section 


stochastic 


geometric 


Experimental 


ratios 


vacuum 


model 


values 


®p-K 1 @pp 


0.69 ±0.02 


0.59 ±0.02 


0.63 


pkI°p-k 


0.83 ±0.08 


0.77 ±0.08 


0.87 



Considering all Hp systems at a given energy, Eqs. ( Jl9| ) and (|20| ) lead to a nonzero 
minimum possible value for the slope, given by 

B^(s) = V a 2 + 1% = l - m 2 + l -B vp {s) = 2.69 GeV" 2 + \b w {s) . (21) 

The existence of this minimum slope that can be observed in the scattering of any hadron by 
a proton is characteristic of the model of the stochastic vacuum. To relate the observables 
for different Hp systems at a given energy, we call G = a[(g 2 FF)a 4 ] 2 a 2 and write 
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B PP - B Hp ( 7 /2)(g p 2 -^) (7/2) a 2 

The last quantity is fixed, for a given energy. As we go from a hadron H to another, we 
obtain in a plot of Bn p against <th p a line from the point representing the observables of the 
pp system to the limit point a = 0, B = B mm given by Eq. (PH). With j3 = 8/3 we have 

r> rjmin 

D _ omin , -°PP -"Hp 1.5 /qo\ 

Br p - B Rp + — a Hp . (23) 

°pp 

Using the data for the pp system at 23.5 GeV we obtain S Hp = 8-59 + 0.014775 a]f . 
This plot is shown in Fig. 2, together with data of the pp, 7rp and Kp systems [[HJ at 
y/s ~ 20 GeV. The limit point is shown inside a square window in the figure. 
In the case of the geometric model we have 

B^(s) = R 2 p = ~B pp (s) (24) 

and a straight line in a plot of Bn p against a" H p 

B pp - -Bh p = 1 / 25 n 

^pp-^Hp K P #) V2 ' 

g is fixed by pure pp data putting <7h p = 0, Bn p = B pp /2, and using Eq. flTBj). Then 

Bu P = \B PP [1 + ^] . (26) 
2 o-pp 

Comparing the two lines shown in Fig. 2 we may tell, for now with some subjective 
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judgement, if either model describes better the data. 



Fig. 2 - Observables for different hadronic systems at 20 GeV. The straight line is the 
prediction of the geometric model, the square window showing the minimum value for the slope 
predicted for small hadrons colliding with protons at this energy. The upper curve and the upper 
square window represent the predictions of the model of the stochastic vacuum. 

The model of the stochastic vacuum predicts that the slope B for the im system at about 
y/s pa 20 GeV is B nn = i]a 2 + 7 2 S% pa 9.6 GeV~ 2 , while the geometrical model predicts 
B n7T = (2/3) r 2 pa 7.5 GeV~ 2 . This is not a trivial difference, as it tests the contribution of 
a nonperturbative QCD effect. 

4. Conclusions 

The results presented in this work exhibit the simplicity of the connections between 
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hadronic high-energy observables determined by the hadronic sizes and stress the importance 
of geometric relations as indicative of properties of nonperturbative QCD dynamics. 

Both the model of the stochastic vacuum and the geometric description of Povh and 
Hiifner give fair account of the present data, although it seems to us that the relations 
provenient of the QCD calculation are more accurate. The differences between the two 
descriptions are important and interesting, and once they can be fully tested by the data, 
may become crucial. The existence or not of the universal term ^r/a 2 , of unique value for all 
energies, representing a pure nonperturbative QCD contribution to hadronic scattering, is a 
question of fundamental importance. Direct hadronic data on hadronic systems with small 
mesons, such as 0p and ipp would be very interesting for the study of nonperturbative QCD 
effects. Hopefully the factory in Frascati will create the opportunity for these studies. 

While the geometric relations of Povh and Hiifner are basicaly empirical, the quantitative 
details (the form of the functional relations and the intrinsic values of parameters) in the 
predictions of the model of the stochastic model determine fundamental QCD quantities 
using only a small amount of data. The model explains the energy dependence of the 
observables in terms of the energy variation of the hadronic sizes, and relates experimental 
quantities for different hadronic systems, exhibiting properties of the extended nature of the 
interaction, which is determined by the structure of the QCD vacuum. 
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FIGURE AND TABLE CAPTIONS 



Fig. 1 - Test of the parameters of Eq. (|l3|), comparing different models, using the experimental 
quantities of the pp and pp systems. For energies up to 62.3 GeV the values of a T are given by the 
parametrization 0"p Om = (21.70 mb)s 0808 and the values of B are those of the pp data. The values 
of i?A and Ca are obtained with the 541 GeV and 1800 GeV data as inputs. The horizontal lines 
represent the constant Ca with the choices A=l, 0.75 (model of the stochastic vacuum) and 0.5 
(geometric model), x 2 represents the average deviation of the five ISR points from the constant 
line. 

Fig. 2 - Observables for different hadronic systems at 20 GeV. The straight line is the 
prediction of the geometric model, the square window showing the minimum value for the slope 
predicted for small hadrons colliding with protons at this energy. The upper curve and the upper 
square window represent the predictions of the model of the stochastic vacuum. 

Table I - Ratios of the pomeron exchange contributions to total cross sections for different 
hadronic systems. The experimental values are taken from the parametrization of Donnachie and 
Landshoff. 



Cross section 


stochastic 


geometric 


Experimental 


ratios 


vacuum 


model 


values 


&pir / @pp 


0.69 ±0.02 


0.59 ±0.02 


0.63 




0.83 ±0.08 


0.77 ±0.08 


0.87 
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